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Abstract 

Determining  the  speed  of  a  ball  as  a  function  of  its  height  after  it  is  launched 
straight  upward  in  the  absence  of  fluid  resistance  is  a  standard  problem  treated 
at  all  levels  of  introductory  physics.  Inclusion  of  drag  in  the  problem  is  seldom 
covered  even  in  the  introductory  course  for  majors.  But,  in  fact,  this  problem 
can  be  solved  via  a  straightforward  application  of  the  work-energy  theorem  for 
either  linear  or  quadratic  drag,  and  instructively  plotted  on  graphs  of  kinetic 
versus  potential  energy. 


1.  Setting  up  the  problem 

Choose  upward  to  be  the  positive  v-axis  and  the  origin  to  coincide  with  the  launch  position  of 
the  ball.  Let  u,  be  the  ball’s  upward  launch  speed,  h  the  maximum  height  that  the  ball  attains 
before  turning  around  and  beginning  its  descent,  and  Uf  the  ball’s  speed  as  it  returns  to  its 
launch  position.  The  problem  consists  in  finding  the  speed  v  of  the  ball  as  a  function  of  its 
height  y  between  launch  and  return. 

In  the  absence  of  drag,  the  solution  can  easily  be  found  by  conservation  of  mechanical 
energy  E, 

K  +  U  =  E  —  constant,  (1) 

where  K  —  \mv2  is  the  ball’s  kinetic  energy  (assuming  the  ball  has  mass  m)  and  U  —  mgy 
is  the  gravitational  potential  energy  (if  the  reference  level  is  taken  at  the  origin).  The  value 
of  the  constant  energy  can  be  conveniently  expressed  either  in  terms  of  the  potential  energy 
at  the  topmost  point,  £/top  =  nigh,  or  of  the  launch  kinetic  energy,  K,  —  \mv2.  Substituting 
these  expressions  into  equation  (1)  and  rearranging  immediately  gives  v(y)  as  desired. 

On  the  other  hand,  when  fluid  resistance  is  present,  a  dissipative  drag  force  FD  acts  on 
the  ball  during  its  motion.  Let  us  initially  consider  the  specific  case  of  a  baseball  moving 
at  ordinary  speeds  of  throwing  in  air  (i.e.,  not  at  the  turbulent  launch  speeds  of  professional 
baseball  pitchers).  In  that  case,  as  further  discussed  in  section  4  below,  the  drag  force  is 
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quadratic  in  speed  [1]  so  that  it  has  a  magnitude  of  the  form  FD  =  bv2,  where  b  is  a  constant 
of  proportionality  (that  depends  on  the  cross-sectional  area  of  the  ball  and  the  density  of  the 
air),  and  a  direction  opposite  to  the  velocity  v  of  the  ball.  If  we  drop  the  ball  from  a  large 
enough  height  that  it  reaches  terminal  speed  l>t,  then  the  gravitational  and  drag  forces  acting 
on  the  ball  must  become  equal  and  opposite,  so  that 

bvj  =  mg  =>•  b  =— f  =>■  Fd  =  mg  ( — ^  .  (2) 

Uf  \vTJ 

(In  general,  buoyancy  can  be  neglected  provided  the  ratio  of  the  fluid  density  to  the  average 
density  of  the  ball  is  much  less  than  unity.  Otherwise  we  can  rigorously  account  for  the 
buoyant  force  by  reducing  the  effective  gravitational  acceleration  g  from  its  vacuum  value; 
specifically,  the  fractional  decrease  in  g  is  equal  to  this  density  ratio  [2].) 

It  is  convenient  to  normalize  the  energies  by  the  kinetic  energy  the  ball  would  have  if  it 
were  moving  at  terminal  speed,  K  \  =  \mvj.  That  is,  we  introduce  the  dimensionless  kinetic 
and  potential  energies,  K.  =  K/Kj  =  (u/u t)2  and  U  =  U/Kj  =  2 gy/vj.  In  terms  of  these 
new  variables,  equation  (1)  for  the  ball’s  motion  in  the  absence  of  air  resistance  becomes 

U  +  K  =  Wtop.  (3) 

When  subject  to  drag,  the  baseball’s  kinetic  energy  after  launch  will  be  smaller  at  a  given  height 
(or  equivalently,  a  given  potential  energy)  than  it  would  be  in  the  absence  of  air  resistance  for 
the  same  launch  speed,  since  the  drag  force  is  dissipative.  We  therefore  need  to  replace  K,  in 
this  equation  by  some  function  of  /C  that  varies  more  slowly  than  linearly.  On  the  other  hand, 
the  function  must  reduce  back  to  /C  in  the  limit  that  /C  <$C  1  because  this  implies  that  v  «  Ur 
which  in  turn  means  that  the  drag  force  is  negligibly  weak  (compared  to  mg)  according  to 
equation  (2).  As  is  formally  proven  in  the  next  section,  the  upward  motion  of  the  ball  is  in  fact 
described  by  replacing  K.  by  the  simple  function  ln(  1  +  1C),  which  by  inspection  satisfies  both 
of  these  conditions.  During  the  downward  phase  of  the  ball’s  motion,  the  drag  force  reverse 
direction  (since  it  is  always  opposite  to  v),  and  the  relevant  function  is  instead  —  ln(l  —  K,). 
(The  insertion  of  two  minus  signs  ensures  that  this  expression  still  has  the  correct  limiting 
value  K,  for  /C  <5jC  1.)  To  summarize,  in  the  presence  of  quadratic  fluid  resistance,  equation  (3) 
becomes 


U  ±  ln(l  ±  K)  =  Utop,  (4) 

where  the  upper  (lower)  signs  refer  to  upward  (downward)  motion  of  the  ball.  Equations  (3) 
and  (4)  are  compared  in  figure  1,  taking  Utop  =  1  in  both  cases.  This  latter  condition  implies 
that  the  baseball  was  launched  at  terminal  speed  in  the  absence  of  air  drag  and  at  V e  —  1  U| 
(i.e.,  31%  faster  than  terminal  speed,  to  allow  for  subsequent  losses)  in  the  presence  of  drag. 
It  returns  to  the  launch  position  with  the  terminal  speed  again,  in  the  absence  of  drag,  but 
with  a  speed  of  only  *J\  —  1  j e\)\  (i.e.,  20%  slower  than  terminal  speed)  in  the  presence  of 
air  resistance.  (See  section  3  for  a  derivation  of  these  numerical  values.)  Experimental  data 
using  video  analysis  of  dropped  balls  [3]  or  balloons  falling  onto  a  motion  sensor  [4]  are  in 
good  agreement  with  the  assumption  of  quadratic  drag  underlying  equation  (4). 

2.  Derivation  of  equation  (4) 

During  an  infinitesimal  change  dv  in  the  ball’s  position,  mechanical  energy  is  lost  due  to  the 
work  done  by  the  nonconservative  drag  force. 


d K  +  dU  =  Fd  •  dr. 


(5) 


Energy-based  solution  for  the  speed  of  a  ball  moving  vertically  with  drag 


1143 


normalized  PE 

Figure  1.  Normalized  kinetic  versus  potential  energies  for  a  ball  launched  straight  up  and 
subsequently  returning  to  the  origin,  assuming  it  is  thrown  just  fast  enough  to  reach  a  maximum 
height  of  Vj/2g.  In  the  absence  of  drag,  the  ball  regains  its  original  speed  on  the  way  down  (dashed 
line),  whereas  quadratic  drag  results  in  a  net  loss  in  kinetic  energy  (solid  curves). 


where  the  ball’s  displacement  is  vertical,  dr  =  dvj.  Noting  that  the  drag  force  is 
Fd  =  ^fjng(v/ur)2j,  where  we  continue  to  use  the  convention  that  the  upper  (lower)  sign 
refers  to  upward  (downward)  motion  of  the  ball,  we  can  rewrite  equation  (5)  as 

dK  +  d U  =  =p ?ig  ^ ^  dy.  (6) 

Despite  the  double  sign  appearing  on  the  right-hand  side,  that  expression  is  always  negative 
because  dy  >  0  when  the  ball  is  travelling  upwards,  while  dy  <  0  when  the  ball  is  descending. 
That  is,  mechanical  energy  is  being  dissipated  in  both  the  upward  and  downward  segments  of 
the  motion,  as  expected. 

As  before,  we  normalize  this  energy  expression  by  dividing  every  term  through  by  Kj  to 
get 

d/C  +  dW  =  TlCdU.  (7) 


which  rearranges  into 
d/C 
1  ±/C 


-dW. 


(8) 


We  integrate  both  sides  over  any  macroscopic  increment  of  motion  (in  either  the  upward  or 
downward  direction)  to  obtain 


±Aln(l  ±/C)  =  -AW,  (9) 

where,  as  usual,  ‘delta’  (A)  refers  to  the  change  in  the  value  of  the  quantity  that  follows  it 
(e.g.,  AW  equals  the  value  of  U  at  the  end  of  the  increment  of  motion  minus  its  value  at  the 
beginning  of  the  increment).  In  particular,  for  the  increment  between  any  arbitrary  point  on 
the  ball’s  path  and  the  topmost  point, 

±  [ln(l  ±  /Ctop)  -  ln(l  ±  K)]  =  -[Wtop  -  If].  (10) 

But  /Ct0p  =  0  since  the  ball  instantaneously  comes  to  rest  at  the  topmost  point  and  therefore 
the  first  logarithm  on  the  left-hand  side  of  this  expression  drops  out.  Rearranging  what  is  left 
now  gives  equation  (4),  as  we  wished  to  show. 


1144 


C  E  Mungan 


3.  Applications  of  equation  (4) 

Apply  equation  (4)  to  the  initial  launch  of  the  ball  upwards  to  find 


ln(l  +/C;)  =  W, 


top 


L/'-p 

h  —  —  In 

2g 


1  +  1  — 
Vj 


(11) 


which  is  the  maximum  height  attained  by  the  ball  in  terms  of  its  launch  speed.  Next  apply 
equation  (4)  to  the  return  of  the  ball  downwards  to  its  launch  position  to  obtain 


-  ln(l  -  /Cf)  =  U, 


top 


l/rr 

h  =  —  In 

2g 


1 


(12) 


(13) 


_1  -  (uf/uT)2 

In  contrast,  equation  (3)  in  the  absence  of  drag  implies  that 

v~  vf 

ZQ  =  Wt0p  =  /Cf  =>■  h  —  —  =  — , 

2  g  2  g 

which  agrees  with  equations  (11)  and  (12)  using  the  Taylor  expansion  ln(l+x)  «  xfoi\x\  <ZC  1 
(as  is  appropriate  because  the  drag  force  can  be  ‘turned  off’  by  taking  the  limit  uj  — >■  oo). 

Equations  (11)  and  (12)  can  be  inverted  to  get  the  normalized  launch  and  return  speeds  in 
terms  of  the  maximum  height. 


Vi 


(2§h 
=  ,  ex P  ~ 

Vr  \  \  V f 


1 


and 


exp 


t  (14) 


Squaring  these  results  and  substituting  h  =  vf/2 g  gives  the  y-intercepts  of  the  solid  curves 
in  figure  1,  consistent  with  the  speeds  mentioned  at  the  end  of  section  1.  In  general, 
dividing  the  second  equality  in  equation  (14)  by  the  first  and  squaring  gives  the  compact 
expression  ZCf/ZC;  =  exp(— (ftop).  On  the  other  hand,  if  one  equates  the  expressions  for  h  in 
equations  (11)  and  (12)  and  rearranges,  one  finds  the  memorable  result  [5] 

1  1  1 

“2  =  +  ~2’  (15) 
vf  vf  vf 

alternatively  expressed  as  l/Kf  =  1  /  K\  + 1  /  K-\  .  Equation  (15)  tells  us  that  Vf  is  always  smaller 
than  both  Vi  and  vj,  that  we  recover  the  drag-free  result  Vf  =  v,  in  the  limit  that  v,  uy, 
and  that  we  return  to  the  ground  at  terminal  speed  if  u,  uy.  Finally,  one  can  calculate  the 
cumulative  amount  of  mechanical  energy  lost  [6]  S  =  E\ost/Kj  in  normalized  form,  as 


£  =  Ki-K-U 


(16) 


since  U\  =  0,  where  the  dimensionless  kinetic  energy  K.  can  be  expressed  in  terms  of  the 
height  y  using  equation  (4).  In  particular,  at  launch  one  finds  /C,  =  exp(Wtop)  —  1  from 
equation  (11).  Note  that  £iost  becomes  thermal  energy  of  the  atmosphere. 


4.  Comparison  with  linear  drag 

For  a  sphere  of  diameter  D  moving  at  speed  v  through  a  fluid  of  density  p  and  viscosity  rp 
the  Reynolds  number  is  defined  as  R  =  p  Dv/ij.  The  drag  force  is  linear  in  v  (Stokes’  law) 
when  R  <  1,  while  the  drag  force  is  quadratic  in  v  (Newton's  law)  when  R  >  1000  [2].  (At 
intermediate  values  of  R ,  one  gets  a  mixture  of  both  forms  that  is  best  treated  numerically. 
Also,  turbulence  arises  when  R  >  2  x  105.)  For  air  at  room  temperature,  p  —  1.21  kg  m-3 
and  i]  —  1.83  x  10~5  kg  m  s-1,  so  that  a  baseball  of  diameter  D  =  7.4cm  is  clearly  in  the 
quadratic  regime  for  speeds  v  greater  than  0.2  m  s~  1  (i.e.,  the  entire  trajectory  except  within 
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normalized  KE  at  launch 


Figure  2.  Normalized  potential  energy  at  the  maximum  height  versus  normalized  kinetic  energy 
at  launch.  The  curves  in  the  presence  of  drag  fall  below  the  straight  line  in  the  absence  of  drag  due 
to  the  dissipation  of  mechanical  energy,  but  asymptotically  approach  the  drag-free  line  in  the  limit 
of  small  launch  speed  (compared  to  the  ball's  terminal  speed). 


2  mm  of  the  top  turning  point).  On  the  other  hand,  if  the  ball  is  a  dust  grain  of  100  /zm  diameter 
falling  at  a  speed  of  0.1  m  s'1,  then  it  is  in  Stokes’  regime.  In  that  case,  we  have  a  drag  force 
of  magnitude  Fq  =  b'v  where  b'  =  3tc  rjD.  At  terminal  speed  up,  the  forces  balance  so 
that  in  analogy  to  equation  (2),  b'vj  —  mg  =>•  F&  —  ;wg(u/ur)  =  mg\pK,.  Therefore, 
equation  (8)  becomes 

dK'  =  -d U.  (80 

1  ±  V)C 

This  is  easily  integrated  by  making  the  substitution  z  =  1  ±  \fiC-  to  obtain 

2A[±V£-ln(l  ±  VlC)]  =  -AW.  (90 

(For  comparison  purposes,  all  equations  numbered  with  a  prime  for  linear  drag  are  analogous 
to  the  corresponding  unprimed  equations  for  quadratic  drag.)  Again  considering  an  increment 
between  any  arbitrary  point  on  the  ball’s  path  and  the  topmost  point,  one  finds  that 

U  ±  2VK.  -  2 ln(l  ±  V/C)  =  Ut op.  (40 

Using  the  expansion  ln(l  +  x)  &  x  —  x1  /2  for  |x|  1,  equation  (40  correctly  reduces  to  (3) 

in  the  absence  of  drag. 

Applying  equation  (40  to  the  initial  launch  of  the  ball  upwards  implies 


2-s/iCi  —  21n(l  +  sJlCi)  —  Ux 


top 


(  Ui  \ 

=>  h  =  — 

—  -  In  1 

i  +  — 

8 

_Vj 

V  bh7  J 

(IT) 


In  the  absence  of  drag,  equation  (3)  implies  that  /C;  =  Wtop;  together  with  the  first  equalities  in 
equations  (11)  and  (IT),  this  drag-free  result  is  plotted  in  figure  2.  The  linear  and  quadratic 
curves  track  closely  with  each  other  up  to  surprisingly  large  launch  speeds;  in  fact,  the  two 
curves  cross  when  u,  =  1.539ut-  Evidently,  a  measurement  of  the  maximum  height  is  not  an 
effective  way  to  distinguish  experimentally  between  linear  and  quadratic  drag  unless  the  ball 
is  launched  with  a  speed  of  at  least  2uj.  For  this  specific  launch  speed,  equation  (11)  predicts 
that  the  ratio  of  the  maximum  height  in  the  presence  of  quadratic  drag  to  the  height  attained 
in  the  absence  of  drag  is  1  In 5  =  0.402,  while  equation  (IT)  predicts  that  the  ratio  of  the 
maximum  height  in  the  presence  of  linear  drag  to  the  height  attained  in  the  absence  of  drag  is 
1  —  5  ln3  =  0.451. 

The  analogue  of  equation  (12)  can  next  be  written  down  for  the  downward  return  of  the 
ball  and  that  resulting  expression  for  li  equated  to  the  second  equality  in  (1 T)  to  obtain 


fj  +  Vf 
Vj 


=  In 


ut  +  Uj 
Ut  —  Uf 


Uf 

—  =  1  +  VT 

Vj 


1  +  Uj/ur 
exp(l  +  Ui/uT). 


(150 
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where  W[x]  is  the  Lambert  W  function  [7],  There  are  two  negative  solutions  of  this  function 
for  values  of  u,/u[  >  0  and  one  should  take  the  solution  with  the  smaller  absolute  value  to 
ensure  that  Vf/vj  >  0.  (The  other  solution  merely  gives  the  launch  velocity  Vf  =  —u,.)  The 
appearance  of  the  Lambert  W  function  nicely  complements  another  recent  work  in  which  this 
function  also  arises  in  connection  with  a  projectile  subject  to  linear  drag  [8]. 
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